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The prevalent Landau damping theory for classical plasmas does not fully explain the Langmuir 
wave decay in dense plasmas. A dielectric function theory adapted from the condensed matter 
physics is extended in order to be applied to dense plasmas and warm dense matter. This theory, 
accounting for the Umklapp process, predicts much higher decay rates than the Landau damping 
theory, which is in better agreement with the existing experimental data obtained from the metals. 
It is demonstrated that this strong plasmon decay leads to the existence of a parameter regime 
where the backward Raman scattering is unstable while the forward Raman scattering is stable. By 
means of the backward Raman compression, intense x-ray pulses can be created in this regime. 

PACS numbers: 52.35.-g, 71.45. Gm, 42.55.Vc, 78.70.Ck 
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I. INTRODUCTION 

Understanding the warm dense matter and dense plas- 
mas has recently become more important in the context 
of various practical problems, including the inertial con- 
finement fusion [l|, 0] • In these media, many interesting 
phenomena observed in classical plasmas would still oc- 
cur, but they get modified due to the quantum diffraction 
and the degeneracy We consider one such phe- 

nomenon, the decay of the Langmuir wave, which plays 
an important role in various applications [7Hl2j|. 

The Landau damping theory is often used to pre- 
dict the plasmon decay in classical plasmas. On the 
other hand, the plasmon decay rate in metals is com- 
puted by complicated quantum theories accounting for 
the local-field correction or the density functional ap- 
proach [H, [3. However, neither of them is sufficient 
to predict the plasmon decay in warm dense matter or 
dense plasmas. In this paper, we introduce a theoretical 
framework to calculate the plasmon decay in dense plas- 
mas by adapting the dielectric function approach used in 
the condensed matter physics. It is suggested that the 
modification of the free-electron wave packet due to the 
presence of ions is considerable, which in turn modifies 
the plasmon decay profiles. It is shown that the decay 
rate of the long wavelength plasmon is much higher than 
the prediction by the Landau damping theory, and that 
the rate is finite and non-negligible in the consideration of 
the Raman compression. It leads us to show that there 
exists a parameter regime where the backward Raman 
scattering (BRS) is unstable while the forward Raman 



scattering (FRS) is stable, enabling the x-ray compres- 
sion without a premature pump depletion. 

The rest of the paper is organized as follows. The ba- 
sic classical and quantum-mechanical Landau damping 
theory is summarized in Sec. II, and why the preva- 
lent Landau damping theory is insufficient for metals and 
dense plasmas is argued in Sec. III. A dielectric function 
formalism, adapted from the condensed matter physics 
research is introduced and extended to be applicable to 
the dense plasmas and warm dense matter (Sees. IV and 
V) The plasmon damping rate is computed in this new 
framework and the existence of a regime in dense plasmas 
where the BRS can be unstable while the FRS is stable 
is demonstrated (Sec. VI), and then it is concluded in 
Sec. VII. 



II. LANDAU DAMPING THEORY 

Consider a wave of the form 0(x, t) = </>exp(ik • x — 
icot) + <fi* exp(— ik • x + icot). The classical analysis of the 
Landau damping rate, 7 c i, for small e<j) is [15| 



7d 



7T< df(u/k) 

2k 2 dv ' 



(1) 



where k (u>) is the wave vector (frequency), 
\J Aim e e 2 / m e is the plasma frequency, n 



pc — 

(m e ) is the 



electron density (mass) , and / is the electron probability 
distribution function satisfying J fd 3 v = 1. The damp- 
ing rate 7 can be alternatively derived in the dielectric 
function formalism, 



* sku@cims.nyu.edu 

t Current Address: 28 Benjamin Rush Ln. Princeton, NJ 08540 



7 



1 



Im [e(k,w)] , 



(2) 



2 



where e is the well-known dielectric function obtained 
through the random phase approximation; 

e(k,u) = 1 + — — / - — — — . (3) 



ui — k • v 



In the regime where the x-ray compression would be fea- 
sible through the BRS, n e is so high that the electron 
quantum diffraction cannot be neglected. In such a case, 
the quantum version of the dielectric function, the so- 
called Lindhard function [l6j], needs to be used; 



e(k,w) 



where Ej 



f 



, , 2 



k 2 h J u-(E f -Ei)/h 
5 v 2 /2 (E f = m e v 2 /2) is 



d 3 v, (4) 



the electron 

kinetic energy of the momentum m e Vi = m e v — Kk/2 
(rn e Vf = m e v + ftk/2). In order to derive Eq. we 
first solve the time-dependent Schroedinger equation us- 
ing a perturbation analysis for each free electron wave 
packet and arrive at 



|q> = |q)o 



e0exp(ik • x — iut) 
' hu - E(q + k) + E(q) 

e<p* exp(— «k • x + iuit) 



|q + k> 



q-k) 



Huj - E(q - k) + E(q) 

where E(p) = fi, 2 |p| 2 /2m e , | )o denotes the free electron 
wave function, and | ) denotes the perturbed wave func- 
tion. The perturbed electron density is given as 

d 3 q 



Sn e (x,t) 



[S(q)l(x|q)| 2 



where g(q), satisfying J g(q)d 3 q/(2ir) 3 = n e , is the oc- 
cupation number. <5n e (k, w) reads 

d 3 q 



Sn e (k, w) = e<j) 



/(q + k)-/(q) 



2tt) 3 3 hu - E{q + k) + E{q) 



The dielectric function, to the first order in (/>, is 

Aire 2 



(5) 



ek,u) 



f 



k 2 



-a(k,w), 



where a(k, w) = 5n e ()t, w) / ecj> is the susceptibility. The 
imaginary part of Eq. @ is 



w m e f E f 



Et — Ei 



)(f(E f )-f(E l ))d 3 v. 



Note that the integration is done in the velocity space 
which is simply related to the momentum space by 
hq/m e = v. In the limit h = 0, the above equation 
is reduced to Eq. (pj. 



(Eq. ([T])) or the degenerate Lindhard dielectric function 
(Eq. 0). The rate in rare dense plasmas is accurately 
predicted by these approaches, However, the rate in dense 
plasmas, warm dense matter, or metals is not. For an in- 
stance, the comparison between experimentally obtained 
damping rates in Al and theoretical computations exhibit 
large deviation (Fig. [TJ . The prediction by the Lindhard 
dielectric function [16|, |l7[ , Eq. (Q} , is essentially zero in 
this range of k, and is not visible in the figure. The 
damping rate from the DuBois' theory [l4j], accounting 
for the dynamical correlations, predicts the damping to 
be proportional to (k/kp) 2 in the long wavelength limit 
where k goes to zero. However, the data from the electron 
stopping experiments in metals [H, [H| suggest that the 
rate is finite in this limit, which is inconsistent with the 
prediction of the aforementioned theories [20I423T ] . The 
experimental decay rate for k < 0.5kF can be estimated 
as 



va(k) = r](k)uj pc / sec . 



(6) 



where r)(k) = r;o + dr\j dk 2 {kjkF) 2 . For a typical metal, 
0.02 <f] < 0.2 and di]/dk = ar] , where 2 < a < 10 
It is tempted to explain that the finite plasmon decay rate 
for k = arises from the electron-ion collisions. However, 
it is difficult to estimate the decay rate by this argument. 
Furthermore, the electron-ion collisions get diminished in 
degenerate plasmas by the electron diffraction and the 
degeneracy [3, H[ , and it is not clear how the quantum 
effect can be incorporated into the plasmon damping. 

The electron wave packets in metals are distorted from 
those of the free electrons through the strong interaction 
between the electrons and the ion lattices. Adler [24| 
computed the electronic dielectric function in the pres- 
ence of the metal ion lattices, by expanding the wave 
packet in terms of the reciprocal lattice wave vector 
and using the quantum random phase approximation. 
Sturm [2Cj and Hasegawa [25| applied this dielectric func- 
tion formalism to Eq. @ and obtained the Langmuir 
wave decay rate, which is in good agreement with the 
experimental data (l8j . In particular, in the long wave- 
length limit, the damping rate is predicted to be finite. 
The dominant mechanism of the decay in Alkali metals is 
shown to be the Umklapp process; an electron inside the 
Fermi shell gets excited by the Langmuir wave outside 
the shell, obtaining the momentum of k + G, where k is 
the Langmuir wave vector and G is the reciprocal vector. 



III. LANGMUIR WAVE DAMPING IN METALS 
AND BREAKDOWN OF THE LANDAU 
DAMPING THEORY 

The Landau damping rate in the free-electron plasma 
can be computed by either the classical dielectric function 



In dense plasmas or warm dense matter, the electron 
density is so high that many physical properties are sim- 
ilar to those in metals, and so is the decay profile of the 
plasmon damping; the ion density is also high, and the 
distortion of the electron wave packet and the Umklapp 
excitation become significant. 
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FIG. 1: Comparison of various damping rates of a plasma of 
7i e = 1.8 x 10 23 cm -3 , which is a reproduction of Fig. 1 in 
Ref. [2f|. The prediction from the classical dielectric func- 
tion (solid line) for T e = 20 eV is compared with that of the 
DuBois' theory (dashed line) for T e = 0. The experimen- 
tal measurements for AL are adap ted from Fig. 3 in Ref. [53] , 
which are originally from Ref. [2Sf (crosses) , Ref. [2?| (circles) , 
and Ref. [3y] (triangles). 



IV. DIELECTRIC FUNCTION IN METALS AND 
ITS EXTENSION TO DENSE PLASMAS 

In this Section, the dielectric function theory developed 
for the metals is briefly reviewed and extended for the 
case of dense plasmas and warm dense matter. In both 
cases, the electron eigenstates are assumed to deviate 
strongly from the free-electron wave packets, \a), where 
a is an index. In the presence of the potential of the form 
<f>(x,t) — 4>exp(i'k ■ x — iu)t) + 4>* exp(— ik • x + iut), the 
wave packet is modified to be 



*> = k> + E 



E 



ftw — E a , + E c 



|(7i}((7i|exp(ik • x)|ct) 



e<f>* 



—ftw + E ai - E c 



|cri)(cri| exp(-ik • x)|cr), 



where the perturbation is assumed to be weak and the 
perturbation theory of the first order is used. | ) denotes 
the original eigenstate and || ) denotes the perturbed 
eigenstate. The perturbation in the density obtained 
from the above relation is 



Sn(k, uj) 



where 

£(01,<7 2 ) 



E^kOIki II exp(ik-x) || en) 

E ^(/Vl) - /(o-2))/3(o-l,CT2,k,w), 



(er 1 |exp(-ik • x)|ct 2 )(ct 2 | exp(«k • x)|cn) 



The dielectric function, up to the first order in e/>, is given 
to be 



e(k,w) = 1 + 



47rn e e 2 
k 2 



*Kf{vi)-ffo)) 



^ (enj exp(-ik ■ x)|g 2 ) kg | exp(ik ■ x)|ai) 
ftw - E ai + E a2 

With an appropriate choice for the eigenstate in a given 
condition, Eq. (171) can be used in various situations. For 
example, Sturrnpol [27| used the following eigenstate 



k> = k) + 



El*) 



(qil^k) 

E q - E qi 



where q and qi are the wave vectors, and (qi|^|q) is the 
pseudo-potential which is experimentally measurable. It 
should be noted that, in case of metals, (qi | V"|q) becomes 
finite only when qi — q is a reciprocal vector. 

Now we apply the above formalism to dense plasmas 
or warm dense matter, where unlike metals, the ions can 
move as individual particles. Since the ions move much 
slower than the electrons and the Langmuir wave fre- 
quency is much faster than the ion relaxation time, we 
assume that the ions are spatially frozen at Xj. We first 
consider the dynamical property of the electrons in the 
presence of the frozen ions, and then obtain the average 
dynamics by averaging over the probability distribution 
of Xj, following the Born-Oppenheimer's approximation. 
Without loss of generality, X^ can be assumed to have the 
following correlation average: 



<E 

i,3 



exp(zs ■ (X, - Xj)) , 
V 



n/(s), 



where V is the volume of the region under consideration, 
and nj(s) is the static two-point correlation function of 
the ions. For independent ions, n.j(s) = ni. Only one 
species of ions of the charge Z is considered. When the 
ion positions are fixed, the electron's free wave eigenfunc- 
tion is modified to be 



k) = |q)+E 



cZ 3 qi exp(-iqi • Xi)J7(|q- qi|) 



(2tt) 3 



E(q) - E( qi ) 



(8) 

where U(q) is the Fourier transform of the ion-electron 
potential. In the absence of the electron screening, 



f7(q) = 4nZe 2 / 
ceptibility 



|q| 2 , and Eqs. ([7J and © give the sus- 



a(k, oj) = a rpa (k, w) + a d onsc(k, to), 
where a r pa, the well-known Lindhard susceptibility, is 



ftiO — E„ 



E„ 



d 3 q /(k + q)-/(q) 
(2tt) 3 ftw - E{k + q) + E(k) ' 



(9) 



4 



and the adcnse, the dense plasma correction, is 



DEGENERATE ELECTRON CASE 



, d 3 s f d 3 q_ U 2 (s 



(2tt) 3 ly ' J (27r) 3 A 2 (q,k,s) 
/(|q + k + s|)-/(q) 
Hu> - E(\q + k + s|) + £(q)' 



(10) 



where A(q, k, s) is 



^ _1 (q,M) 



1 



£(q)-£(q + s) 
1 



E{q + k) - £7(q + k + s) ' 



(11) 



The imaginary part of ayensc is given as 



Im [a de 



rf 3 s 



r n/(s) 



rf 3 c 



t/ 2 (s) 



(2tt) 3 Jw y (2tt) 3 A 2 (q,k,s) 
xtt(5(Sw - E(q + k + s) + E(q)) 
x(/(q + k + s)-/(q)). 



(12) 



Some comments are made on Eq. (|10[) . which is the main 
result of this paper. First, the integration over s corre- 
sponds to the summation of the contribution from the 
Umklapp process in metals. An electron with an initial 
momentum hk absorbs the energy huj and the momentum 
H(k + s). This absorption is accompanied by an exchange 
of the momenta between an electron and an ion. Second, 
in metals, the integration over s is replaced by the sum- 
mation over the reciprocal lattice vectors, which repre- 
sents the summation of the inter-band transitions in the 
case of metals. Third, the ion correlation is incorporated 
into nj (s). The ions tend to be strongly correlated in the 
warm dense matter, which could be easily accounted for 
by this formalism. Fourth, from Eq. @, the decay of the 
Langmuir wave is given as 



1 1 47re 2 

-Im [e] = - -p-Im (a rpa + a dc 



7 



(13) 



The first term in the right hand side accounts for the 
Landau damping, and the second term is the correction 
due to the presence of the ions. Lastly, o^ense does not 
arise purely from the quantum-mechanical effects, as this 
term does not vanish in the limit h = 0. In the limit 
h = and k = 0, the decay rate becomes proportional 
to nj and T e 3 ^ 2 , which is reminiscent of the electron- 
ion collision frequency. However, the exact coefficient, 
which is derived in Eq. (|12p. cannot be obtained from 
the ion-electron collision argument. When h 7^ 0, the 
quantum degeneracy and the electron diffraction effects 
are self-consistently incorporated into Eq. (TT2"1) . 



Consider the dense plasma where the electrons are 
completely degenerate. The Fermi-energy is given as 
E F = 36.4 x (n/n 24 ) 2/3 eV, where n 24 = 10 24 cm~ 3 . The 
computation in this Section is valid as long as T e < Ep. 
In the partially degenerate and classical case, the electron 
screening and the quantum diffraction need to be treated 
carefully, which would be considered elsewhere [3l[ . 

The susceptibility a™ of the free electron plasma, 
computed by Lindhard [16], is given as 



a rpa (k,w) 



3n e 
m e v F 



h(z, u), 



where vf = \j2Ep /m e is the Fermi velocity, z = k/2kp, 
u = uj/kvp, and h — h r + ihi. The real part of h, h ri is 



h r = 



[l-(z- u) 2 ) log 



II 



82 

-i(l-(z + u ) 2 )log 
v V \z + u + 1| 



z — u + 1 



and the imaginary part hi is 



hi = —u, for |z + u| < 1, 



82 



(l-(z-u) 2 ), for \z - u\ < 1 < \z + u\, 



= 0, if \z-u\> 1, 



In the limit k = 0, due to the delta function in Eq. (|12[) . 
A _1 (q, k, s) in Eq. (fTTjl can be approximated as 



A- 1 (q,k,s) = 



ks 



using the screening potential U(s) = 47rZe 2 /(s 2 + k 2 f ), 



where ktf = V3uj pc /vf is the Thomas- Fermi screening 
length. In this case, A is independent of q, which can be 
taken out of the q-integration in Eq. (fTUj) . Using Eq. ([9]), 
Eq. (0 reads 



Im [a dense (k,w)] 



2 ,2 



d 3 s n 4 (k-s) 2 (47rZe 
(2^) 3 m 2 (fiw) 4 (s 2 +k 2 f ) 2 

xn/(s)Im [a rpa (k + s, w)] , 



which can be further simplified to 

AE 2 Z(k F e 2 ) 2 {An) 2 



Im [a d cnsc (k,w)] = 



(hu^ 3tt 2 J 4(2tt) 3 
(k • s) 2 

Im [a rpa (k + s,u)] , 



k 2 (s 2 + k 2 f ) 2 



(14) 



5 



0.07 




FIG. 2: The damping rate of a degenerate hydrogen plasma 
for a range of the electron density, 10 24 cm~ 3 < n e < 
10 26 cm- 3 . 
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FIG. 3: The damping rates of a degenerate hydrogen plasma 
as a function of the wave vector for two cases of electron 
density. 



where we assumed n/(s) = n/. iy(0) in Eq. ([6]) from the 
above equation is shown in Fig. [2] 

For non-zero fc's, Eq. (ITi|) cannot be used, as A is de- 
pendent on q. A rather complicated expression for A 



A' 1 (q,k,s) 



h 2 k-s 



m e (fajj) 2 
fik- (q + s/2) 



moj 



fik- (q-s/2) 



mjjj 



which should remain in the integrand for the q- 
intcgration in Eq. (|10|). For given s and k, we integrate 
over q first and then over s. For simplicity, we assume 
k = (k x , 0, 0) and s = (s cos 9, s sin 0,0), and denote q = 
qixi + q 2 x 2 +93X3, where i\ — (k + s cos 6*, s sin#, 0)/|k + 
s\, x 2 — (— s sin#, 1 -I- cos9, 0)/|k + s|, and 2 3 = (0, 0, 1). 
In this coordinate system, A depends only on k, s, 9, and 

Q2- 



A- 1 Ck,s,q 2 ) 



h 2 k-s 

m e (huj) 2 



;K K 



where 



sk cos 9 



m e uj 



sk sin 9q 2 sk cos ( 



The real part of a remains still complicated, but the 
imaginary part given in Eq. (IT2"j) could be simplified as 
the delta function would eliminate q\ integration and the 
integrand of dq is independent of q%. We do not present 
the detailed steps here. After some tedious manipula- 
tion, the q-integration is reduced to an one-dimensional 



integration 

d 3 q U 2 {s) 
(2tt)3 A 2 (q,k,s) 



Im 



/(q + k + s)-/(q) 



hw - E(q + k + s) + E(q) 



(ft 2 k-s) 2 f/ 2 (s) fV 1 -^-") 



m 2 (hio) 4 



47T 2 ft 2 |s + k| 



xk + k _1 (R~(kp, z, u, q 2 ) — R + (kp, z, u, q%j) , 



where R ± (kp, z,u,q 2 ) 



^k 2 F {l~(z±u) 2 )-q 2 if 



k 2 F (l - (z ± u) 2 > 0, R ± (k F ,z,u,q 2 ) = if k 2 F (l 
(z ± u) 2 < 0, z = |k + s|/2fejr, and u = Lu/kvF- The 
right-hand side of the above equation is a function only 
of k, s, and 9. We do the s-integration in the spherical 
coordinate and arrive at 

2irs 2 dnds (h 2 ksn) 2 U 2 (s) 



Im[a(k,o;)J = J 



-nj(s) 



m 2 (huj) 4 



kf^y l — (z—u) 2 



dqi 



47r 2 ft 2 |s + k 



x (R~{k F ,z,u,q 2 ) - R + {k F ,z,u,q 2 )) 



This is a three-dimensional integration that needs to be 
evaluated numerically. In Fig. ([3]), r\ in Eq. ([5]) as a func- 
tion of k is shown for a hydrogen plasma of two different 
electron densities. 



VI. RAMAN SCATTERING 

In this Section, we discuss the existence of a regime 
where an x-ray pulse can be compressed by the BRS. 
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More detailed presentation can be found in Ref. [26| for 
the case of metals. Here we provide a similar computa- 
tion in the case of degenerate dense hydrogen plasmas. 

The BRS is a three-wave interaction where two light 
pulses and one Langmuir wave exchange the energy I32i— 
34 1 . A one-dimensional interaction can be written as [33l . 
3 = 



Vl 



v 2 



d_ 

dx 
d_ 
dx 
d_ 

dx 



V2^ M = 

v-i ) A 3 = 



-iciA 2 A 3 , 
-ic 2 AiAl, 
-ic 3 A 1 A 2 ', 



(15) 



where A; = eEi/m e LOiC is the ratio of the electron quiver 
velocity of the pump pulse (i = 1) and seed pulse (i = 2), 
relative to the velocity of the light c, and A 3 — n e /n e is 
the intensity of the Langmuir wave. i>i (v 2 ) is the rate 
of the inverse bremsstrahlung of the pump (seed), v 3 is 
the plasmon decay rate, Cj = uj 2 e /2uJi for i = 1 and 2, 
C3 = (cq) 2 /2u> 3 , and LUi(ui2) is the frequency of the pump 
(seed). In the BRS (FRS), the energy conservation is 
given as wi = cj 2 + u pe (uji = uj 2 + w pe ) and the momen- 
tum conservation is given as q = k\ — k 2 (q = k\ + k 2 ). 

The inverse bremsstrahlung (v\ and v 2 ) in metals is 
given as [H, |36[ , 



v{w, E) = AirniZ 2 



F(a) 



ra 2 v% w 2 



2 K i/2' 



(16) 



where k = (ftiu /2m e v F ) < 1, Dp is the Fermi energy and 
2 - 2e 2 E 2 /m 3 Hu} 3 with E being the electric field of the 



a 



pulse, and the function F(a) can be approximated as 
a 2 /6 when a < 2 and (2/ira) log(2a) 2 otherwise [35[. In 
Eq. (|16j) . it is assumed that Hlj > Ef and fno > T e where 
T e (Ef) is the electron temperature (Fermi energy). The 
above equation can be simplified, when a < 2, as 



2.9 x 10 15 Zk 



i , 2 

-3/2 ^0 



/ sec . 



(17) 



If 1/v in Eq. (fT7| is shorter than the pulse duration, the 
pump pulse will be heavily damped before being com- 
pressed. Note that v\ is proportional to oj 2 c /ui 2 . 

From Eqs. ^ and (fTS")) . the BRS growth rate, assum- 
ing that the pump intensity is large enough to make the 
BRS unstable, is roughly estimated to be 



A\ uj 2 



Upe/ sec. 



(18) 



pc 



The larger gs is, the stronger the BRS is. 

Now, a possible operating regime for the Raman com- 
pression is obtained from Eqs. (|6|), (flB"]) . and (fT7|). The 
stability condition for the FRS is given, from Eq. flT51) , as 
C2C3I | 2 < ^2^3, where v 2 is obtained from Eq. (|17p and 



v 3 is from Eq. ©, and c 3 — uj pc /2. This condition can 
be written as |^4i| > |^4if| where 



\A 



IF 



^V 2 (qF)v 3 bJ 2 



(19) 



and qp = cj pc /c. The same condition for the BRS can be 
estimated, using c 3 = (cji + uj 2 ) 2 /2uj pc , as \Ai\ < |^4ib|, 
where 



I , 1 2 v 2 (q B )v z 4w pc w 2 
\A\b\ - 2 7 i 72 



(20) 



and q B = (wi + ui 2 )/c. \A 1B \ 2 in Eq. (|20J) is 
smaller than |^4if| 2 in Eq. (|19| by a factor of 
(wi + w 2 ) 2 /w 2 e (?7(g F )/77(g B )), which is very large since 
vi ( lF)/i r l{ ( lB) — 1- Therefore there exists a regime where 
the BRS is unstable but the FRS is stable. This is due 
to the strong inverse bremsstrahlung and a rather slowly 
varying damping of the Langmuir waves as a function of 
k. 

We provide an estimation for hydrogen plasmas with 
n e — 10 24 cm~ 3 and 10 26 cm -3 . In the former case, when 
hujp C = 37.11 eV and uJi/ujp C = 10, it is estimated that 
r](k) = 0.06 from Fig. [3] The threshold for the FRS esti- 
mated from Eq. (fT9j) is |Aif| 2 = 3x 10 -4 , which corre- 
sponds to 3.4 x 10 19 W/cm 2 . The threshold for the BRS 
estimated from Eq. is |^is| 2 = 0.0075xl0~ 4 or8.4x 
10 16 W/cm 2 . For n e = 10 26 cm^ 3 , hujp C = 371.19 eV, 
and wi/wp = 10, 77(0) = 0.18 and r](k = 2ui/c) S 0.25 
from Fig. [3] The threshold for the FRS estimated from 
Eq. (HI]) is |^if| 2 = 0.288 x 10~ 4 which corresponds 
to 3.225 x 10 20 W/cm 2 . The threshold for the BRS es- 



timated from Eq. ^DJ) is |Ai B | 2 
1.116 x 10 18 W/cm 2 . 



0.000998 x 10~ 4 or 



VII. CONCLUSION 

The prevalent Landau damping theory predicts the 
plasmon decay rate to vanish in the long wavelength 
limit. On the other hand, the observations from the elec- 
tron stopping experiments in metals suggest that it is 
finite in that limit. Using a well-known dielectric func- 
tion formalism developed in the condense matter, we 
propose a theoretical framework predicting the plasmon 
decay rate in dense plasmas more accurately (Eq. (fTTJ)) ) 
than conventional theories. A few advantages of our ap- 
proach are in order: First, the robustness of this frame- 
work is verified through the experiments with metals [2fj| . 
Second, the quantum degeneracy and the diffraction of 
the electrons are self-consistently accounted for. Third, 
the inter-ion correlations could be incorporated for a 
wide range of physical parameter regimes through nj(s), 
the static two-point correlation function among the ions. 
Lastly, the theory takes into account of the Umklapp pro- 
cess. The computation based on this dielectric function 
approach suggests that the plasmon decay rate in dense 
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plasmas in the long wavelength limit is much higher than 
the prediction by other theories. 

This strong decay has implications for various pro- 
cesses in dense plasmas, warm dense matter and metals. 
For instance, it renders the FRS less harmful for the Ra- 
man compressor in metals and dense plasmas, compared 
to the FRS in rare dense plasmas. The plasmon decay 
in metals strongly depends on the direction relative to 
the lattice structure. As the electron temperature T e in- 
creases, the electron screening of the ions gets weaker, 
resulting in a stronger plasmon damping [20f, H3] ■ As T e 
increases further, the ions would finally lose their lattice 
structure. A larger value of 7]{k) might be preferred for 
a strong pump, as it would weaken the FRS; however, a 
small value of r)(k) might be preferred for a weak pump, 
as a strong BRS is necessary. 

As a consequence, there exist a parameter regime in 
metals where the BRS compression is strong while the 
FRS is stable (Sec. VI). As an example, a pump pulse 
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